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ABSTRACT 


This  paper  presents  a  statistical  software  package 
developed  for  use  on  the  Apple  II  Plus  aicrocoiputer, 
modified  with  the  Apple  Pascal  language  card.  The  program 
addresses  the  following:  determination  of  confidence 
intervals  for  single  and  bivariate  populations;  hypothesis 
testing  for  one  and  two  paraaetars;  coaputation  of 
cuaulative  distribution  values  for  the  Normal,  Student's  T, 
Chi-square#  P#  Binoaial  and  Poisson  distributions; 
coaputation  of  quantile  values  for  the  Noraal,  Student's  T, 
Chi-square  and  F  distributions.  The  program  also  has  the 
capability  to  store,  retrieve,  and  modify  data  for  use  with 
the  statistical  procedures.  The  program  was  written  in  (FCSD 
Pascal,  which  because  of  its  portability  indicates  that 
little  or  no  modification  would  be  required  to  use  it  with 
other  computers  which  are  BCSD  Pascal  compatible.  In 
addition,  because  of  Pascal's  block  structure,  the  program 
can  be  easily  modified  or  enhanced  to  include  other 
statistical  procedures  which  are  of  interest  to  the  user. 
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I.  INTRODUCTION 

For  the  operations  analyst  or  for  that  Batter  anyone  who 
utilizes  the  aethodology  and  problem  solving  approach  of  the 
operations  analyst,  soae  form  of  computational  device  is  a 
necessity.  C03t,  rigid  interface  requirements,  and  a  number 
of  other  factors  have  in  the  past  frustrated  and  stifled  the 
analyst  in  bringing  to  bear  the  necessary  coiputational 
power  on  his  problem.  In  the  late  1970's,  however,  Texas 
Instruments  introduced  a  significant  amount  of  computing 
power  packaged  in  the  programmable  TI-59  calculator. 
Revolutionary  is  perhaps  a  bit  too  strong  as  a  description 
of  the  impact  that  this  and  like  devices  have  made  in  the 
operations  research  community,  but  certainly  most  analysts 
would  agree  that  the  amount  of  computing  power  that  can  now 
be  held  in  one's  hand  and  used  to  solve  problems  is  indeed 
remarkable. 

In  spite  of  this  impact,  however,  the  hand-held 
calculator's  contribution  might  even  now  be  waning  and 
yielding  to  a  more  spectacular  capability  found  in  the 
microcomputer.  The  growth  in  the  capability  of  these 
devices  since  the  beginning  of  the  industry  in  1971  has  been 
phenomenal.  Nastrakas  details  this  growth  and  gives  a 
glimpse  of  the  possible  direction  that  the  1980 's  will  see 
in  this  industry  [Ref.  ij  .  One  factor  accounting  for  the 
growth  in  microcomputer  technology  and  one  that  serves  to 
insure  its  future  is  the  intense  competition  that  pervades 
the  industry.  When  one  thinks  of  hand-held  calculators,  he 
thinks  of  Texas  Instruments  or  perhaps  Hewlett-Packard.  When 
one  thinks  of  microcomputers,  he  may  think  of  hpple,  Exidy, 
North  Star,  PET,  TIS-30,  or  a  number  of  other  devices  with 
like  capability. 
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A  typical  aicrocoaputer  package  costing  between  $2,000 
and  $2,500  eight  consist  of  the  following:  the  computer 
(about  the  size  of  a  seall  portable  typewriter)  with  64  K 
bytes  of  usable  randoe  access  memory,  two  floppy-disk  drives 
(5  1/4  inch  dimeter)  for  prograe  storage,  and  a  black  and 
white  eonitor  for  output  display.  4  variety  of  prograeaing 
languages  are  also  available  including  the  University  of 
California,  San  Diego  (OCSD)  Pascal,  BASIC,  FORTRAN ,  PILOT, 
etc.  soee  of  the  nicrocoaputers  also  have  a  graphics 
capability  allowing  the  user  to  visualize  aatheaatical 
foras,  plot  graphs,  and  plot  data  observations. 

The  Aaerican  populace  is  being  conditioned  through 
current  periodicals  and  news  features  to  expect  an 
increasingly  iaportant  role  for  the  aicrocoaputer  in 
everyday  life  —  froa  grocery  shopping  to  environaental 
control  for  the  hoae.  T he  Deceaber  1,  1980  issue  of  "o.  S. 
News  and  Norld  Report'*  predicts  that  eighty  percent  of  the 
United  States*  households  will  have  a  aicrocoaputer  by  1990. 
In  view  of  the  prospective  proliferation  of  these  devices 
and  also  the  significant  computing  power  that  now  exists, 
the  operations  research  analyst  can  ill  afford  not  to  begin 
to  exploit  the  capabilities  of  aicrocoaputers.  Indeed,  the 
hardware  capabilities  have  grown  and  are  growing  so  rapidly 
that  today  good  compatible  software  to  support  these 
impressive  capabilities  is  seriously  lacking  in  specialized 
fields  such  as  operations  research.  It  is  this  software 
deficiency  that  this  thesis  addresses. 
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II.  BACKGROU  ND 

The  operations  analyst  can,  on  occasion,  be  required  to 
establish  confidence  intervals  or  test  hypotheses  about  an 
unknown  parameter  from  a  known  or  assumed  population.  The 
computer  program  written  as  a  part  of  this  thesis  allows  the 
analyst  to  quickly  and  easily  accomplish  these  tasks  when 
observations  are  from  normal.  Exponential  or  Bernoulli 
populations. 

A.  DESCRIPTION  OP  THE  MICROCOMPO TE B  SYSTEM 

The  software  development  was  lone  on  an  Apple  II  Plus 
microcomputer  with  two  floppy-disk  drives  (5  1/4  inch 
diameter)  as  add-on  peripherals.  The  system  is  equipped 
with  a  language  systea  giving  the  capability  to  use  the 
University  of  California,  San  Diego  (UCSD)  Pascal 
programming  language  a3  well  as  the  BASIC  language  which 
comes  resident  with  the  computer.  The  standard  Apple 
computer  has  a  forty  column  output  display  which  makes  it 
compatible  for  use  with  a  standard  television  set.  This 
capability  can  be  enhanced  to  an  eighty  column  display  with 
an  additional  peripheral  device,  provided  a  monitor  is  used 
for  display  in  lieu  of  a  television.  The  output  format  for 
the  program  is  written  for  an  eighty  column  display  device; 
however,  using  the  special  built-in  features  of  the  Apple, 
the  format  can  easily  oe  made  to  display  split  screen  on  a 
systea  not  equipped  with  an  expanded  display  peripheral. 

3.  SOFTWARE  DEVELOPMENT 

Ling  and  Muller  give  several  considerations  which  should 
be  observed  in  the  development  of  software  for  statistical 
analysis  [Ref.  2,  Ref.  3].  Among  these  considerations  are 
the  following: 
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1.  Choice  of  Programing  Language 

The  Pascal  programing  language  used  in  this 
software  development  effort  has  a  auober  of  features  which 
make  it  particularly  attractive  for  use  with  microcomputers. 
First#  Pascal  is  a  very  concise  language.  The  compiler  is 
saall  and  compact  and  fits  easily  within  the  available  off¬ 
line  storage  space  of  the  floppy-disk.  Second,  Pascal  is  a 
high-level,  general-purpose  language  |Ref.  '•*]  .  The 
language  was  originally  introduced  in  1971,  which  is  recent 
in  comparison  to  aost  high-level  languages.  It  was  intended 
to  be  used  in  teaching  new  prograaaers  good  techniques  and 
style.  The  Pascal  language  fully  exploits  the  fundamental 
concepts  of  structured  programing,  which  is  a  technique 
used  by  many  professional  prograaaers  to  write  large  complex 
computer  programs  |Ref.  hj  .  The  use  of  these  techniques 
facilitates  developing  programs  in  a  modular  fashion  (i.  e. , 
break  the  overall  pactage  into  logical  sub-packages  and 
proceed  to  program,  debug,  and  validate  each  sub-package 
individually) .  The  final  step  in  the  process  is  to  combine 
the  sub-packages  to  form  the  overall  package.  Osina  the 
attendant  statistical  package  as  an  example,  six  sub¬ 
packages  or  modules  comprise  the  overall  package.  Each 
module  is  independent  of  the  others  and  can  stand  alone  when 
compiled  with  the  main  program.  Third,  the  Pascal  language 
performs  arithmetic  computations  significantly  faster  than 
the  BASIC  language.  Pascal  is  often  implemented  as  a 
“pseudo"  interpreted  language  meaning  that  the  text  versions 
of  programs  are  first  compiled  into  a  code  file.  It  is 
during  the  compilation  phase  that  syntax  errors  are  detected 
by  the  compiler  and  brought  to  the  programmer's  attention. 
This  code  file  is  interpreted  and  executed  during  the 
execution  phase  of  the  program.  Rost  processors  can  and 
typically  do  interpret  a  Pascal  code  file  significantly 
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faster  than  a  corresponding  BASH  progras  which  performs  the 
saae  computations  ^Ref .  4|  . 

2.  Computational  Efficiency 

In  spite  of  thair  impressive  capabilities, 
microcomputers  are  decidedly  infarior  to  largar  computers  in 
the  two  key  araas  of  computational  speed  and  accuracy. 
Pascal,  as  implemented  on  the  Apple,  will  only  perform 
computations  using  six  decimal  places  of  accuracy  ^Eef.  5] 
and  displays  only  five  places  past  the  decimal.  Accordingly, 
algorithms  which  are  used  and  work  very  well  on  larger 
computers  might  have  no  chance  of  producing  the  same  results 
on  microcomputers  simply  because  the  execution  time  is 
excessive,  or  they  require  double  precision  arithmetic. 


C.  PORTABILirr 

A  portable  program  is  one  that  can  be  run  on  a  number  of 
different  computer  systems  ^Ref.  2]  .  Programs  written  in 
the  Pascal  language  are  valuable  in  this  regard  in  that  they 
may  be  run  on  a  variety  of  microcomputers  without 
alteration.  Some  of  the  machines  for  which  this  is  true  are 
computers  with  the  following  microprocessors:  3080,  3085, 
Z80,  6502(Apple),  6800,  and  9900.  Portability  is  a  very 
strong  asset  of  the  Pascal  language  even  though  it  is 
accomplished  at  the  expense  of  reduced  computational  speed. 
The  compiled  code  version  of  programs  is  called  “ p-code". 
Bach  machine  has  a  special  interpreter  program  which  takes 

li  «  . 

the  p-code  and  converts  it  to  a  form  compatible  with  the 
existing  host  xicroprocessor . 

D.  BASIL  I  RECRI37A3LB  I SFORM ATI0S  PILES 

Ling  suggests  the  use  of  help  files  to  allow  the  user  to 
make  efficient  use  of  a  program  ^Ref .  2] .  The  number  of 
information  files  which  are  available  and  their  content  is 
clearly  a  matter  of  judgment  based  on  assumptions  concerning 
the  knowledge  of  the  user  population  about  the  procedures 


used  ia  the  program.  Anticipating  that  most  of  the  users  of 
the  accompanying  statistical  package  will  be  familiar  with 
the  basic  concepts  of  the  procedures  themselves,  only  one 
help  file  which  pertains  to  data  entry  requirements  is 
included. 


3.  SISPL2  USES  INTERFACE 

The  user  is  called  upon  to  make  numerical  entries 
throughout  the  prograa.  Because  Pascal  is  a  strongly  typed 
language,  a  variable  of  type  integer  cannot  be  assigned 
floating  point  values.  Indeed,  this  requireaeat  is  so 
strict  that  if  the  prograa  expects  the  user  to  enter  an 
integer  and  he  accidently  enters  a  number  in  decimal  point 
notation,  the  program  will  abort  and  cause  the  entire  system 
to  re-initializ e.  situations  such  as  this  are,  of  course, 
undesirable.  Soreover,  all  programs  written  with  the  intent 
of  establishing  a  dialogue  between  the  prograa  and  the  user 
should  be  as  trouble-free  as  possible  for  the  user  and 
ainiaize  as  much  as  possible  the  user's  chances  of 
committing  a  fatal  error  when  responding  to  program  prompts 
or  entering  data.  In  general,  schemes  to  accomplish  this 
are  costly  in  terms  of  computational  efficiency  and 
programming  steps;  however,  user  convenience  is  almost 
always  worth  the  costs.  The  problem  with  data  entry  or 
numerical  entry  from  the  keyboard  is  addressed  in  this 
statistical  package  by  making  all  numerical  entries  using 
a  string  and  converting  the  string  to  the  numerical  value 
it  represents.  A  string  is  a  variable  type  which  is  a 
linear  array  of  characters.  For  example,  given  a  string 
variable  called  'S',  which  is  assigned  a  value  of  '235', 
the  character  in  array  position  *  s[l]'  is  '2',  'S^2  '  is 
•3',  and  'S^j'  is  *5'.  Since  the  characters  'O'  through 
'9'  have  corresponding  numerical.  American  Standard  Code 
for  Information  Intercaange  (ASCII)  values  of  forty-eight 


through  fifty-seven,  conversion  is  accomplished  by 
subtracting  forty-eight  froa  the  character's  ASCII  value, 

A  Pascal  procedure  which  converts  strings  to  umbers  is 
shown  in  Appendix  B. 

The  advantages  in  using  a  schaie  such  as  this  for  all 
numerical  entries  are  two-fold.  First,  the  user  can  take 
advantage  of  the  direct  cursor  addressing  available  on  many 
microcomputer  video  displays  to  correct  a  data  entry  prior 
to  entering  it  into  coiputer  memory.  This  is  not  possible 
when  the  program  expects  real  numbers  as  input.  Second,  the 
program  segment  which  converts  the  string  may  return  either 
an  integer  or  a  real  number,  whichever  is  required  by  the 
program.  This  avoids  having  the  user  concerned  with  the 
typing  requirements  demanded  by  the  Pascal  language. 


Ill 


Algorithms  used  for  statistical  computations  on 
microcomputers  should  be  selected  with  the  goal  of  providing 
the  best  accuracy  achievable  at  the  minimum  computation 
time.  Given  the  limitation  of  only  six  decimal  places  of 
accuracy  on  the  Apple  microcomputer#  many  algorithms 
requiring  more  precision  in  their  computations  must  be 
rejected.  The  algorithms  remaining  as  candidates  must  be 
carefully  screened  to  insure  that  their  required 
computations  are  reasonable  from  a  time  standpoint  and  that 
they  exercise  fully  the  accuracy  capability  of  the  machine. 


A.  CALCULATION  OP  VASIA  MCE 

To  illustrate  the  complications  stemming  from  reduced 
accuracy,  consider  the  following  example  given  by  Forsythe 
[Ref .  6]  . 


Pind  the  variance  of  the  following  set  of  numbers: 
48499,  48503,  48500,  48498,  48500 
The  definition  for  the  variance  is 


S  * 
where 

This 
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»  -  1 


(X,  -  ?)* 


S*  =  variance 

f=  number  of  observations 
=  sample  mean. 


formula  can  be  expanded  to  the 
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N  -  1 


following  form: 


Use  of  this  formula  in  the  Apple,  however,  would  yield  a 
variance  of  zero,  simply  because  the  required  accuracy  is 
not  available  for  computations.  The  answer  of  zero  is,  of 
course,  incorrect.  Forsythe  in  his  article  on  statistical 
computing  offers  the  following  alternative  algorithm  to 
compute  the  variance. 


1  3 


sua  :*  0.0; 

S2  ;=  0.0; 

READ  (N)  ; 

POR  I  :s  1  TO  S  DO 
BSGIH 

READ  (X)  ; 

DEVIATION  :=  X  -  SOM ; 
son  :=  SUM  ♦  DEVIATION/! ; 

S2  :=  S2  ♦  DEVIATION  *  (X  -  S(JH|  ; 

END; 

S2  ;=  S2/(N  -  1) ; 

This  algorithm  produces  the  correct  answer,  "variance  * 
3. 5**.  It  also  illustrates  how  many  of  the  limitations  of 
the  microcomputer  can  be  overcome  through  careful  selection 
of  algorithms. 

B.  DISTRIBUTIONS  AND  INVERSES 

The  algorithms  used  to  compute  probability  distributions 
and  inverses  and  the  source  of  each  are  listed  in  Appendix 
A.  Nhea  compared  to  the  tabular  values  listed  in  Dixon  and 
Nassey  |fief .  7]  ,  the  algorithms  are  accurate  to  at  least 

three  decimal  places  ia  probability  with  the  exception  of 
the  P  distribution.  The  P  distribution  is  accurate  to  three 
decimal  places  in  probability  in  almost  all  cases;  however, 
some  values  may  differ  from  the  listed  tabular  values  by  as 
much  as  .002  in  probability.  Although  the  ?  quantiles  may 
differ  slightly  from  the  listed  tabular  values,  the 
probabilities  corresponding  to  the  values  given  by  the 
algorithm  are  accurate  to  three  decimal  places.  All  of  the 
algorithms  produce  results  within  three  seconds  except  the  T 
distribution.  Nhen  computing  foe  very  large  degrees  of 
freedom  for  the  T  distribution,  computation  time  is  a 
function  of  the  degrees  of  freedom.  Typically  1000  degrees 
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of  freedom  takes  approximately  eight  seconds. 

Often,  increased  accuracy  can  only  be  gained  through 
iterative  techniques.  This  is  the  case  with  many  of  the 
probability  distributions  and  their  inverses,  the 
T  distribution  being  an  exaiple.  In  the  algorithms  shown  in 
Appendix  A  much  of  the  excessive  computation  time  has  been 
alleviated  by  combining  two  or  more  algorithms.  Por 
example,  in  the  Chi-square  distribution,  the  P  distribution 
and  the  inverse  F  distribution,  one  algorithm  is  used  for 
small  degrees  of  freedom  and  another  for  large  degrees  of 
freedom.  Good  algorithms  for  these  distributions  exist  for 
large  degrees  of  freedom  which  are  not  based  on  iterative 
techniques  and  hence  are  computationally  fast.  However,  for 
small  degrees  of  freedom,  their  accuracy  falls  off  rapidly. 
Conversely,  the  algorithms  using  iterative  techniques  are 
very  accurate  at  all  ranges,  but  3low  for  the  larger  degrees 
of  freedom  since  the  number  of  iterations  required  is 
proportional  to  degrees  of  freedom.  The  break  between  small 
and  large  degrees  of  freedom  is  purely  subjective,  based  on 
choosing  the  best  combination  of  speed  and  accuracy. 
Selection  of  the  algorithms  themselves  was  likewise  based  on 
the  best  combination  of  speed  and  accuracy. 

C.  POISSON  AND  BIN3MIAL  DISTR IBUTIDNS 

Large  mainframes  have  the  capability  of  computing  the 
Poisson  and  Binomial  Distributions  directly  from  their 
definitions. 


POISSON: 

If  the  random  variable  I  is  distributed  Poisson  with 
parameter 

ji?  j! 


then  <=  kj 


avx 
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BINOMIAL: 


If  the  randoi  friable  X  is  distributed  Binomial  with 
parameters  {n,pf, 

then  p[x  <=  k]  =  |*j  p*  (I-?)"'3 

Both  of  these  expressions  contain  factorials  and 
summations;  consequently  for  large  values  of  'k’  in  the 
Poisson  and  1  n'  and  in  the  3iaomial,  execution  time 

light  be  excessive,  or  the  intermediate  values  in 
computation  light  exceed  *10  •  (the  maximum  nuaber  capable 

of  being  represented  oa  the  Apple).  A  better  solution  for 
the  calculation  of  these  probabilites  is  to  use  the 
Chi-square  identity  for  the  Poisson  and  the  P  identity  for 
the  Binomial  £aef .  8,  Paf.  9]  . 

POISSON: 

Given  X  ~  Poisson 

then  Pr  [x  <*  k]  *  1  -  X1(xk*i) 

where  2k  ♦  2  are  the  degrees  of  freedom 
of  the  Chi-square  variate. 

BINOMIAL : 

Given  X  ^  Binomial  (n,pj 

t“~  *  ir  <*  kl  *  f  <=  H*-  • 

where  2n-2k  equals  the  degrees  of  freedom  for  the  numerator 
and  2k*2  equals  the  degrees  of  freedom  for  the  denominator. 

When  k  *  n,  the  probaoility  is  1.0. 

D.  STATISTICAL  ALGORITHMS 

Derivation  of  the  algorithms  used  in  finding  confidence 
intervals  and  hypothesis  testing  is  shown  in  Storer 
[aef .  10] . 
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The  options  available  to  the  user  are  shown  in  the  block 
diagram  (Figure  1)  .  ilhen  the  prograa  is  executed,  the  user 
begins  in  the  “outer  level".  To  proceed  he  mst  select  froa 
one  of  che  seven  options  which  ace  shown  on  a  “menu"  (Figure 
2)  . 

A.  THE  DATA  ENTRY  30D0LE 

Selecting  option  '5'  froa  the  outer  level  aenu  (Figure 
2)  will  cause  the  data  entry  aodule  aenu  (Figure  3)  to 
appear  on  the  screen. 
i •*  general 

The  data  entry  requirements  are  intended  to  be  as 
trouble-free  as  possible  for  the  user.  The  user  is  prompted 
for  data  input  by  the  following  line: 

Record  N  — > 


•  N'  is  an  integer  sequentially  updated  by  the 
program  when  the  'return1  key  is  pressed.  Poliowing  the 
arrow,  the  user  inputs  as  aany  data  values  as  he  wishes 
with  entries  separated  by  one  or  aore  spaces.  The  only 
restriction  is  that  he  should  not  exceed  the  length  of  the 
display  line.  The  nomenclature,  record,  indicates  only  a 
logical  or  convenient  grouping  of  data  from  the  user's 
point  of  view. 

Prior  to  entering  any  data,  the  user  is  asked 
whether  or  not  the  observations  he  is  entering  are  paired. 
Since  two  of  the  statistical  procedures  are  predicated  on 
paired  data,  answering  'yes'  to  this  question  will  cause 
suaaary  statistics  to  be  coaputed  on  'X,I'  pairs.  These 
suaaary  statistics  are  only  good  for  use  with  the 
procedures  requiring  paired  observations.  All  data  entered 
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figure  1.  bloeK  biugrain 


WtKHIiM. 


1)  Hypothesis  Testing  One  Parameter 

2)  Hypothesis  Testing  Two  Pacaneters 

3)  Confidence  Intervals  Single  Population 

4)  Confidence  Intervals  Bivariate  Populations 

5)  Distributions  and  Inverses 

6)  Data  Entry 
Q)  uit 

Figure  2.  The  Outer  Level  Menu 


1)  instructions 

2)  Create  a  new  data  file 

3)  Correct/Add  to  existing  datafile 

4)  Enter  data  without  storing 

5)  Review  existing  datafile 

6)  Review  suaaary  statistics  of  existing  file 
Q)  uit  and  return  to  outer  level 

Figure  3.  The  Data  Entry  Hodule  Menu 
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last  be  in  ‘I,!'  pairs  with  the  'x*  observation  listed 
first.  All  other  statistical  procedures  and  their  attendant 
data  entry  requirements  assume  that  all  of  the  observations 
come  from  a  single  population.  These  general  requirements 
are  available  to  the  user  vhen  he  selects  option  '1'  from 
the  data  entry  menu  (Pigure  3) . 

2.  Create  A  Hew  Data  Pile 

Selecting  option  *2*  (Figure  3)  will  prompt  the  user 
to  specify  a  fill  name  for  the  data  observations.  After 
responding  ml*.*  ?.  file  name,  the  operating  system 
establishes  *  **>;«ctory  entry  at  the  beginning  of  the 
largest  anwfcfcl  £lock  of  space  containing  at  least  fifteen 
blocks  the  specified  disk.  Because  the  filing  system  in 
the  University  of  California,  San  Diego  (UCSD)  Pascal 
language  random  access,  each  data  file  entered  will  have 
allocated  fifteen  blocks  of  space  to  insure  that  there  is 
enough  room  to  extend  the  file  if  the  user  desires  to  do  so 
at  some  later  time. 

The  floppy-disk  (5  1/4  inch  diameter)  used  by  the 
Apple  System  provides  a  storage  space  of  280  blocks.  This 
results  in  a  capability  to  store  seventeen  data  files  on 
each  disk.  Each  data  file  can  contain  a  maximum  of  ninety 
records.  Since  each  record  can  contain  as  many  observations 
as  the  user  desires  provided  that  it  does  not  exceed  the 
length  of  the  display  line,  a  reasonable  planning  figure  is 
eight  data  observations  per  record.  This  results  in  an 
upper  limit  of  approximately  720  data  entries  per  data 
file. 

Data  entry  for  a  new  file  begins  with  ’Record  1 1 . 

The  user's  only  concerns  vhen  entering  data  should  be  to 
separata  each  observation  by  one  or  more  spaces,  to  not 
exceed  the  length  of  the  display  line,  and  to  enter  'X,X' 
pairs  if  he  has  previously  indicated  paired  observations. 
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Since  the  data  is  entered  as  a  string  variable  and  then 
converted,  it  is  possible  to  correct  any  entry  prior  to 
going  to  the  next  record  by  simply  using  the  *  back  arrow' 
key  to  move  to  the  place  at  which  the  correction  is  to  he 
sale.  The  Legal  syabols  which  nay  be  used  in  data  entry  are 
tea  digits  'O'  through  '9',  the  decimal  point,  the  comma, 
the  plus  (+|  and  minus  {-)  signs,  and  the  '  E'  for  scientific 
notation.  Commas  are  ignored  by  the  procedure  which 
converts  the  strings  and  are  included  only  as  a  convenience 
for  the  user. 

Pressing  the  'return'  key  at  the  and  of  a  record 
teninates  that  record  and  prompts  the  user  to  input  the 
next  record  in  sequence.  If  the  user  inadvertantly  enters 
an  illegal  character  while  entering  data,  the  program  will 
advise  hii  of  this,  indicate  what  the  character  was,  and 
proipt  him  to  reenter  the  record. 

Once  the  user  has  entered  all  of  his  data,  he  must 
press  the  'escape'  key  and  the  'return'  key  following  tne 
last  data  entry.  This  will  close  the  data  file.  when  the 
file  is  closed,  the  following  summary  statistics  are 
computed  on  the  observations  in  the  file: 

a.  The  Sample  dean 

b.  The  Sample  Standard  Deviation  (N-1) 

c.  The  Sum  of  the  Observations 

d.  The  Sum  of  Squares  of  the  Observations 

e.  The  Number  of  Observations 

For  paired  data,  the  same  statistics  are  computed; 
however,  they  are  computed  on  the  differences  of  the 
pairs.  Hence,  the  number  of  observations  for  a  data  file  of 
paired  observations  is  exactly  half  of  the  total  number  of 
data  entries. 
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The  suaaary  statistics  for  each  data  file  are  kept 
in  a  separate  file  that  is  initially  established  on  the  same 
disk  as  newly  created  data  files.  The  file  of  summary 
statistics  requires  one  block  of  disk  storage  space.  It  is 
differentiated  froa  the  original  file  of  observations  by  an 
'S'  concatenated  to  the  original  file  naae. 

3.  Correct/idd  to  Existing  Data  File 

Selecting  option  '3*  (Figure  3)  while  in  the  data 
entry  aodule  will  proapt  the  user  to  specify  a  file  naae. 
Once  the  data  file  is  retrieved,  the  aonitor  will  show  the 
naae  of  the  file,  the  auaber  of  records  in  the  file,  and  the 
nuaber  of  observations.  Iaaediately  after  this  information 
is  the  following  proapt: 

Enter  Record  Humber  — > 

Selecting  any  record  number  between  '1'  and  the 
total  nuaber  of  records  in  the  file  will  cause  the  retrieval 
of  that  record  and  will  display  as  follows: 

Old  Record: 

Record  2  — >  12345 

You  aay  replace  the  complete  record  or 
press  <RTM>  for  no  changes. 

Record  2  — > 


Pressing  'return'  leaves  the  existing  record 
unchanged  and  prompts  the  user  to  enter  another  record 
nuaber.  If  corrections  are  to  be  lade,  the  user  must  '.nter 
the  cooplete  new  record  opposite  the  lower  proapt  and  press 
'return'.  The  updated  record  will  display  as  follows: 

Record  2  — >  12335 

Press  <RT N>  if  OK,  <250  if  not. 

Pressing  'return'  completes  the  update  and  prompts 
for  a  new  record.  Pressing  'escape'  produces  the  following 
display: 
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Enter  corrected  data  aad  <HTN> 
Record  N  — > 


The  user  then  may  retype  the  line.  This  sequence 
■ay  he  repeated  until  the  record  appears  as  the  user  wants. 
Entire  cecocds  are  erased  or  deleted  by  typing  a  space  to 
produce  a  blank  line.  If  the  user  desires  only  to  update  an 
existing  data  file  and  not  extend  it,  he  types  the  nuabar, 

' -1',  in  response  to  tha  proapt  for  record  number,  whictx 
closes  the  updated  file. 

If  the  user  selects  a  record  nuaber  that  is  greater 
than  the  nuaber  of  records  in  the  file,  he  reenters  the  data 
entry  phase  beginning  with  the  record  nuaber  iaaediately 
following  the  last  existing  record.  Por  example,  if  there 
ace  thirty-seven  existing  records  in  the  file,  selecting  any 
nuaoer  greater  than  thirty-seven  will  produce  the  following 
display: 

Enter  new  records. 

Enter  <ESC>  as  the  last  entry  and  <RTM >  to  terminate  input. 
Record  38  — > 

From  this  point,  the  user  proceeds  exactly  as  if  he 
were  in  the  data  entry  phase  and  terminates  by  pressing  the 
•escape*  key  iaaediately  following  the  last  entry. 

4.  Enter  Data  Without  Storing 

This  option  allows  the  user  to  compute  the  mean, 
standard  deviation,  sua  of  observations,  sub  of  squares  of 
observations,  and  nuaber  of  observations.  The  format  for 
entering  data  is  exactly  the  saae  as  previously  discussed, 
for  creating  a  new  file.  is  indicated,  data  is  not  stored 
under  this  option;  hence,  once  a  record  is  terminated,  the 
entries  for  that  record  cannot  be  recovered  nor  changed. 

This  option  provides  an  expedient  way  to  determine  the 
suaaary  statistics  of  a  group  of  observations. 


5.  Review  an  Existing  Data  File 

Selecting  this  option  alloys  the  user  to  quickly 

review  any  or  all  of  an  existing  data  file  in  blocks  of  ten 
records  at  a  tine.  The  user  is  prompted  to  enter  the  file 
naae  of  an  existing  file.  The  program  retrieves  that  file 
and  displays  the  first  ten  records.  pressing  'return*  at 
this  point  causes  the  next  ten  records  to  be  displayed, 
etc.  Pressing  the  'escape*  key  at  the  end  of  any  display, 
returns  the  user  to  the  menu  for  the  data  entry  module 
(Figure  3) . 

6.  Review  summary  Statistics  of  an  Existing  File 

All  files  that  are  stored  on  disk  have  an  associated 
summary  statistics  file  that  is  created  by  the  program  when 
a  newly  created  or  updated  data  file  is  closed.  This  file 
contains  the  following  information: 

a.  dean 

b.  Standard  Deviation 

c.  Sum  of  Observations 

d.  Sum  of  Squares 

e.  Number  of  Observations 

The  summary  statistics  fila  only  is  called  when 
specifying  a  data  file  to  be  used  in  the  other  modules 
containing  the  statistical  procedures.  Because  of  this,  it 
is  not  necessary  to  have  the  original  data  file  on-line  when 
performing  the  statistical  procedures;  only  the  summary 
statistics  file  is  required.  When  using  the  filing  system 
resident  in  the  Apple  Pascal  language  system  to  obtain 
directory  listings  of  various  disks,  the  sumaary  statistics 
file  is  distinguished  by  a  concatenated  'S'  on  the  ^nd  of 
the  original  file  name.  For  example,  if  the  original  file 
name  was  STAT:0ATA1,  then  the  corresponding  summary 
statistics  file  is  named  STATtDATAlS.  Since  the  length  of 
all  summary  statistics  files  is  one  block,  it  is  possible  to 
access  a  maximum  of  27 ’A  files  of  summary  statistics  on  any 
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one  disk  of  storage. 

7.  Q)  ait 

Pressing  the  ' Q*  key  will  return  the  user  to  the 
outer  level  (Figure  2). 

B.  DISTRIBUTIONS  AND  20ANTILES  MODULE 

Selection  of  any  of  the  options  froa  the  distributions 
and  guantiles  aodule  aenu  (Figure  4)  will  produce  further 
proapts  which  require  the  user  to  enter  the  necessary 
inforoation  concerning  values  of  the  randoa  variable, 
degrees  of  freedom,  and  probabilities  as  appropriate. 

After  each  computation  the  following  prompt  appears: 

C)  ontinue  or  Q)  uit 

Pressing  the  'C'  key  will  allow  the  programmer  to 
continue  calculation  in  the  previously  selected  distribution 
or  quantile.  Pressing  the  *Q*  key  will  return  the  user  to 
the  aenu  for  this  aodule  (Figure  4) . 

The  Pascal  language  system  allows  the  user  to  develop 
his  own  specialized  libraries  of  often  used  subroutines  for 
general-purpose  or  special-purpose  computations.  It  is  in 
such  a  library  that  the  algorithms  for  the  distributions  and 
quantiles  are  kept.  Using  a  special  library  has  two  major 
advantages.  First,  when  the  user  is  developing  the  aaia 
program,  he  is  not  penalized  by  extra  compilation  time  for 
any  of  the  routines  in  the  library.  The  code  in  the  liorary 
is  linked  by  a  separate  process  ^Ref.  hj .  Second,  the 
liorary  can  be  used  by  other  programs  which  require  the  use 
of  the  algocithas  contained  therein.  The  algorithms  for  the 
distributions  and  quaatiles  fit  logically  in  a  library  since 
it  is  likaly  that  other  statistical  packages  will  require 
taair  use.  The  Apple  reference  manual  explains  the 
procedure  used  to  establish  new  libraries  j^Sef.  1 1  j  . 
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1 

2 

3 

4 

5 
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Normal  Distribution 
T  Distribution 
Chi-Square  Distribution 
F  Distribution 
3inoaial  Distribution 
Poisson  Distribution 


A) 

B 

C 

D 


Normal  Quantiles 
T  Quantiles 
Chi-Square  Quantiles 
F  Quantiles 


Q)  uit  and  return  to  outer  level 


Pigure  4.  The  Distributions  and  Quantiles  Nodule  Menu 


POPULATION  ASSUMPTIONS 


1 

i 

5 

6 
Q 


Population 


Assumptions 


Noraal 
Normal 
Normal 
Normal 
Bernoulli 
Exponential 

uit  and  return  to  outer  level 


u?. sigaa-sqr  known 
u  S  sigaa-sqr  ? 
u  &  sigma-sqr  ? 
u  known,  sigaa-sqr  ? 


u  3=>  mean 

sigma-sq  ==>  variance 


Parameter 

u 

u 

sigma-sq 

sigma-sq 

mtIf 


Pigure  5.  Confidence  Intervals  Single  Population  Menu 
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C.  CONFIDENCE  INTERVALS  AND  HYPOTHESIS  TESTIN3  SODULES 

Selecting  either  the  confidence  intervals  nr  hypothesis 
testing  options  will  cause  the  lanu  for  that  soiule  to  be 
displayed.  All  of  these  lenus  contain  infonation  siailar 
to  that  shown  in  figure  5  for  confidence  intervals. 

The  parameters  about  which  an  interval  is  to  be  computed 
or  a  Hypothesis  is  to  oe  made  are  listed  in  the  right  hand 
column  of  the  menu.  The  distribution  from  which  the 
observations  came  and  assumptions  about  the  populations  are 
listed  in  columns  one  and  two,  respectively. 

1 .  Data  Requirements 

When  one  of  the  options  is  selected  that  does  not 
involve  Bernoulli  or  Poisson  populations  the  following 
display  appears: 


1 

2 

3 

4 


Use  existing  data  file 
Enter  data  and  store  it 
Enter  data  w/o  storing 
Use  summary  statistics 


Becausa  of  the  nature  of  the  observations,  any  tests 
or  intervals  involving  Bernoulli  or  Poisson  observations  are 
entered  using  summary  statistics;  hence,  for  these  cases, 
this  display  is  skipped. 

a.  Use  existing  Data  File 

The  user  is  prompted  for  the  name  of  the  data 
file.  When  the  file  name  is  entered,  the  program  retrieves 
the  summary  statistics  file  associated  with  that  file  from 
disk.  In  the  case  of  bivariate  populations,  two  file  names 
are  needed,  the  first  containing  the  • X '  observation  and  the 
second  containing  the  'I'  observation.  Paired  observations, 
as  noted  previously,  are  entered  in  one  file. 

b.  Enter  Data  and  Store  it 

Selecting  this  optioa  will  display  a  short 
message  to  the  user  informing  him  that  all  data  storage  must 
be  accomplished  from  the  data  entry  module  (Figure  3) .  The 
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user  has  the  option  at  this  point  to  enter  the  data  without 
storing  it  or  to  return  to  the  outer  level  (Figure  2)  and 
select  the  data  entry  module. 

c.  Enter  Data  a/0  Storing  It 

Data  is  entered  in  the  same  format  as 
previously  discussed  in  the  section  pertaining  to  the  data 
entry  module.  No  permanent  disk  record  is  made  of  the 
entries.  Hence,  once  each  record  of  observations  is 
terminated,  there  is  no  way  to  retrieve  it  to  make 
corrections.  Where  bivariate  populations  are  required,  the 
user  is  prompted  to  enter  all  of  the  'X*  observations  first 
and  all  of  the  '7*  observations  second.  Por  paired 
observations  all  of  the  'X,Y'  pairs  are  entered  as  one 
populati on. 

d.  Use  Summary  Statistics 

On  all  othar  tests  or  intervals,  use  of  summary 
statistics  is  optional  except  as  previously  mentioned  for 
data  from  Bernoulli  or  Poisson  populations.  &t  each 
prompt,  the  user  is  asked  only  for  the  necessary 
information  to  perform  the  statistical  procedure  he  has 
selected.  The  distinction  is  made  in  each  prompt  whether 
or  not  the  statistics  required  are  the  sample  parameter 
values  (sample  mean,  sample  standard  deviation)  or  the  true 
parameter  values  (true  mean,  true  standard  deviation) . 

2.  Confidence  Intervals 

When  computing  confidence  intervals,  the  user  must 
supply  the  additional  information  concerning  the  desired 
level  of  confidence  and  the  type  of  interval  (two-sided, 
one-sided  upper,  one-sided  lower).  Typically,  computation 
of  the  desired  interval  takas  one  or  two  seconds  and  is 
displayed  as  follows: 
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95.00  Percent  confidence  intervals  for  u 

Saaple  dean  3  3.000 

Standard  deviation  3  1.581 

Upper  a  4.964 
Lover  3  1.036 

Another  interval  using  sane  data,  Y)es  M)  o  — > 

Following  the  display,  the  user  is  asked  whether  or 
not  he  desires  to  compute  another  interval  using  the  saae 
data.  If  he  responds  'yes',  he  aay  then  vary  the 
confidence  level  and/or  the  type  of  interval  without 
having  to  again  specify  the  data  base.  Answering  ‘no*  will 
return  the  user  to  the  menu  for  the  module  in  which 
computations  are  currently  being  performed. 

3.  Hypothesis  Testing 

Hypothesis  tests  require  the  user  to  specify  the 
null  hypothesis.  Typically,  the  hypotheses  involve  'equal 
to',  'less  thaa  or  equal  to',  or  'greater  than  or  equal  to' 
comparisons  and  are  displayed  for  the  user  in  a  fora 
similar  to  the  one  below: 

1)  u  3  u  [0l 

2)  u  <=  u  CO] 

3)  u  >3  u  [0] 

The  symbol  [o]  represents  the  null  hypothesis 
value.  The  user  enters  this  value  if  required  by  the  test. 
An  example  display  following  computation  is  as  follows: 

HYPOTHESIS:  U  3  U  [01 

Saaple  mean  3  3.033 

u  COT  3  2.500 

The  P-value  is  0.519 

Another  test  using  the  saae  data,  Y)es  N)  o  — > 


The  user  is  not  told  to  accept  or  reject  the 
hypothesis;  rather,  he  is  given  a  p-value  as  shown  above. 
The  p-value,  or  probability  level,  is  an  indication  of  the 
level  of  confidence  associated  with  the  hypothesis  ^Sef. 
12|.  High  p-values  convey  a  high  confidence  in  the  null 


hypothesis;  conversely,  low  p-values  reflect  a  lack  of 
confidence  in  the  validity  of  the  hypothesis.  &  p-value  of 
.05,  for  example,  indicates  that  if  the  hypothesis  is  indeed 
true,  there  is  only  one  chance  in  twenty  that  the  data  used 
in  the  test  is  consistent  with  the  hypothesis.  Upon 
completion  of  the  test,  the  user  has  the  option  to  perform 
another  test  with  the  same  data  or  return  to  the  menu  for 
the  module  in  which  tests  are  currently  being  performed. 
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The  decade  of  the  1930’s  promises  to  be  particularly 
bright  in  teras  of  affording  the  operations  research  analyst 
easy  access  to  computing  power.  Microcomputers  now 
available  and  their  aore  capable  descendants  will  doubtless 
play  an  iaportant  role  in  securing  this  access.  However, 
equally  as  iaportant  are  the  software  packages  which  will 
accoapany  these  computers. 

The  software  package  described  in  this  paper  provides 
the  analyst  with  a  useful  sat  of  statistical  tools  which  can 
be  used  on  one  of  the  aost  popular,  current  aicrocoaputers, 
the  Apple.  The  five  aajor  modules  contained  in  the  package 
(confidence  intervals  for  single  and  bivariate  populations, 
hypothesis  testing  for  one  and  two  parameters,  and  the 
distributions  and  quantiles)  are  designed  to  be  easy  for  the 
analyst  to  use  and  to  cushion,  as  much  as  possible, 
potential  user  mistakes.  The  algorithms  used  throughout  the 
program  were  chosen  on  the  basis  of  being  compatible  with 
the  microcomputer  with  respect  to  size  and  computing 
precision  and  providing  the  best  combination  of  speed  and 
accuracy. 

Pascal,  the  programming  language  used,  offers  not  only 
the  advantages  of  portability  and  increased  computational 
speed,  but  also  flexibility.  Pascal  is  flexible  in  that 
large  complex  programs  are  programmed  in  modular  segments 
which  are  then  combined  into  the  overall  program.  It 
follows  that  programs  developed  in  this  way  are  easily 
enhanced  by  the  addition  of  new  modules.  Such  is  the  case 
with  this  statistical  package  which  could  be  significantly 
enhance!  by  the  addition  of  a  regression  and  ao  analysis 
of  variance  module. 


This  package  and  those  to  follow  which  ace  compatible 
with  current  and  future  micr oconputer s  can  have  a 
significant  impact  in  the  analyst  community  in  two  key 
areas.  First,  the  analyst  can  be  better  educated.  Simply 
alleviating  the  tedium  which  accompanies  the  application  of 
many  statistical  procedures  will  give  fledgling  analysts  the 
opportunity  to  work  more  problems  and  be  exposed  to  a 
greater  variety  of  situations  in  the  school  environment. 
Perhaps  of  equal  importance,  the  educational  environment  can 
provide  the  opportunity  to  accustom  the  analyst  to  the 
capabilities  that  can  and  should  be  available  for  his  use  in 
a  working  environment.  Second,  by  expanding  computing  power 
into  areas  which  were  not  privileged  to  have  it  before,  the 
educational  process  has  a  better  opportunity  to  continue. 
Today,  it  is  reasonable  to  assume  that  the  professional 
growth  of  many  analysts  is  stifled  from  a  lack  of  computing 
machinery  with  which  to  attack  his  problems. 

Taking  full  advantage  of  the  microcomputer's  hardware 
capabilities  requires  efficient  compatible  software.  In 
specialized  areas  such  as  operations  research,  the  analysts 
themselves  must  logically  provide  the  bulk  of  the  effort  in 
software  development.  The  statistical  package  which  is  the 
subject  of  this  thesis  effort  scarcely  begins  to  provide  the 
full  complement  of  tools  which  the  analyst  requires.  If  the 
operations  research  community  is  to  take  advantage  in  a 
timely  manner  of  the  new  opportunities  afforded  by 
microcomputers,  effort  must  continue  now  in  software 
development. 


32 


APPENDIX  A 


FUNCTION  Z  ; 

!********************************** ***** 
* 

* 

*  NORMAL  DISTRIBUTION 

c*  SOME  COMMON  BASIC  PROGRAMS 

*  3RD  ED.,  P.  128 

'* 

*************************************** 


CONST 
Cl  = 
C2  * 
C3  = 
C4  = 
C5  = 


0. 4361836: 
-0. 1201674; 
0. 937298; 

0.  33267: 
2.50662§3; 


BEAL; 


0.0  THEN 


VAR 

XX,XT,R,T  I 
3EGIN 

IP  STDEV  <* 

BEGIN 
ERROR; 

EXIT  (Z)  ; 

END 
ELSE 
BEGIN 
XT  :  =  X 

x  : = 

xx  :=  x*  . 

R  :=  EXP<-XX/2.  0) /C5; 

X  :=  1.0/(1  .0  ♦  C4*ABS(XJ )  ; 

T  :=  0.5  -a*(C1*X  ♦C2*X*X  +C3*X*X*X) 
IP  IT  <  MEAN  THEN 
Z  : =  0.5  -  T 
ELSE 

Z  :  *  T  ♦  0.5; 

END; 

END;  (*  Z  *) 


ABS  ( (X-MEAN)  /STDEV) 
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FUNCTION  INVZ 


****  ********************  *********  -  ***** 
* 

* 

*  NORMAL  0IJANTIL2S 

*  HANDBOOK  OF  MATHEMATICAL  PONCTIONS 

*  P.  933 

* 

* 

*************************************** 


CONST 

Cl  =  2.515517; 
C2  =  0.302853; 
C3  »  0.010323; 
D1  *  1.432788; 
D2  *  0. 189269 ; 
D3  »  0.  001308; 

VAB 

PT, T, NO  M, DEN  : 


BEAL; 

(P  <=  0.0)  THEN 


BEGIN 

IF  (P  >*  1.0)  OR 
BEGIN 
EBBOB; 

EXIT  (IN7Z)  ; 

END 

ELSE  IF  P  >  0.5  THEN 
PT  :  =  1.0  -  P 
ELSE 

T  :=  SQBT’(LHM.0/(PT*PT)))  ; 

MOM  j a  Cl  ♦  C2*T  ♦  C3*T*T; 

DEN  l 3  1.0  ♦  D1*T  ♦  02*T*T  *  D3*T*T*T; 
IF  P  >  0.5  THEN 

INVZ  :=  T  -  NOH/DEN 
SLS  E 

INVZ  ;=  -  (T  -  NOH/DEN)  ; 

)z  I*  INVZ  *) 


END;  (■ 
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PO NOTION  T 


***  **********  ******  *******  ***«.  > 


T  DISTRIBOTION 
CACN 

ALGORITHM  344 


***  ************** ********************* 

CONST  Cl  =  0.63661977; 

VAR 

I.N  :  INTEGER; 

ANS,D2,P1,F2,T1,T2,XT  :  REAL; 


BEGIN 

IP  DP  <  1 
BEGIN 
ERROR; 
EXIT  (T) 
END 
ELSE 
BEGIN 

XT  ;=  X: 
X  :=  ABS 

T 1  :=  x/ 

T2  : =  1. 


THEN 


=  ABS(X); 

;=  X/SQRT (DP 


T2  :  =  1.0?(1.0+M*P1)  ; 

IP  ODD  (DP)  THEN 
BEGIN 

ANS  :=  1. 0-Cl*ATAN (T1) t 
IF  DP  <>  1  THEN 
BEGIN 

D2  i=  C1*T1*T2: 

ANS  i=  ANS  -  D2; 

END; 

IF  DP  <>  3  THEN 
BEGIN 

FI  :=  0.0: 

N  :=  (OF- 2)  D£V  2; 

FOR  I  :=  1  TO  N  DO 
BEGIN 

P2  ;=  2. 0*1- PI ; 

D2  :*  D2*T2*P2/(P2*1)  ; 
ANS  ;*  ANS-02; 

END; 

END; 

END 

ELSE 

BEGIN 

D2  :=  T1*SQRT  (T2)  ; 

ANS  ;*  1.0-D2: 

IF  D2  <>  2  THEN 
BEGIN 

FI  :*  1.0: 

N  ;*  (DP-  2)  DIV  2: 

FOR  I  ;*  1  TO  N  DO 
BEGIN 

F2  :*  2 . 0 *1  - P 1 ; 

02  :*  D2*T2*F2/  (F2  +  1) ; 
ANS  :*  ANS-D2; 

END;  (*  END  POR  *) 

END; 

END:  l*  END  ELSE  *) 

IP  ANS  <  0.  0  THEN 


:  =  C1*T1*T2: 
i  i=  ANS  -  D2 ; 


THEN 


0.0  : 
(DF-5) 
;=  1 


D£  V 
TO  N 


:*  2.0*I-P1 ; 

:*  D2*T2*P2/(P2*1)  ; 
;*  ANS-02; 


ASS  :=  0.0: 

IP  XT  <  0.0  THEN 
T  :=  ANS/2.0 
ELSE 

T  :*  1.0  -  ASS/2.0; 

SSD; 

END;  (*  T  *) 


FONCTION  INVT 


T  QUANTILES 

C  ACM 

ALGORITHM  396 

******************  ********************* 

* 
* 
* 
* 
* 
* 
* 

*************************************** 


CONST 

HALFPI  =  1.570796327; 

TAB 

PT,DEN,A,B,C,D,X, Y  :  REAL; 

BEGIN 

IF  (DP  <  1)  OB  (P  >=  1.0)  OH  (P  <=  0.0)  THEN 
BEGIN 
EBBOB; 

EXIT  (INVT)  ; 

END 

ELSE  IF  P  >  0.5  THEN 
P  :=  2.0*  (1.0  -  P) 

ELSE 

P  :=  2. 0*P; 

IF  DF  =  2  THEN 

INVT  :=  SQBT(2.0/(P*  (2.0  -  P>)-2.0) 

ELSE 

BEGIN 

IF  DF  =  1  THEN 
BEGIN 

?  1=  P*HALFPI; 

INVT  :=  COS  (PJ/3IN  (P) ; 

END 

ELSE 

BEGIN 

A  :=  1.0/  (DP  -  0.5); 


C  :=  ( ( 20 /do*  A/'B  -  98.0)  *A  -  16.0)*A  +  95.36  : 

D  ;=  ({S'l-S/tB  ♦  C)  -  3 . 0)  /B  ♦  1 . 0)  *S3aT  ( A*HALF?I)  *D? ; 
X  :=  D*P; 


Y  I*  XPN ( X,  (2.0/DP 
IF  Y  >  0.05  ♦  A  THE 
BEGIN 

X  :=  IN7Z  (P*0.  5)  ; 
I  :=  X*X: 

IP  DF  <  5  THEN 


!■  C 


((  (0.05 


*  0.  3*  (DP  -  4.  5)  *(X 
0 5*D*X  -  5.0)  *X  -  7.0)  * 
“  4*Y  *  6.  3)  *Y  +  36.0)  * 


C 

c  :  = 

Y  :  = 

3  +  1.0 

Y  ;* 

IP  Y 

Y 

ELSE 

Y  :*  0. 5*  (Y*Y)  ♦  Y; 

END 
ELSE 

Y 


+  0.6) 

7.0‘ 


X  -  2. 0)  *X  ♦  B  ♦  C; 
Y  ♦  94.5)  /C-Y-3.0)  / 


A*  (i*Y)  ; 

>0.002  THEN 
!■  EXP(Y)  -  1.0 


IP  ?T  >* 
INVT  : 


(( 1. 0/ (( (DP  +  6.0)  /  (DP* 
(DP  +  2. 0)  *3.0)  +0. 5/(D? 
(DF+1.0j/(DF  +  2.0)M.O 


0.5  THEN 

:  SQBT  (DP*Y) 


DF*Yi-0. 089*D-0.  322)  * 
?  . 0)  )  *Y- 1.0)* 

0/Y; 
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ELSE 

INYT  :*  -  SQfiT (D?*Y)  ; 
END; 

2ND ; 

END;  (*  INVT  *) 
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FUNCTION  CHI5Q  ; 

r**«*  **********  ************************** 


*  * 

*  *’ 

*  CHI-SQUARE  DISTRIBUTION  *' 

*  SONS  COMMON  BASIC  PROGRAMS  * 

*  IDF  <*  40)  P.  130  * 

*  HANDBOOK  OP  MATHEMATICAL  FUNCTIONS  * 

*  (DF  >  40)  P.  941  * 

*  * 


(*«**  ************  ***** *******  ******  ****** 
?  AR 

I  :  INTEGER: 

Y, POHER, TEMP, NU3, DEN, J,L,M  :  REAL; 


BEGIN 

IP  (DF  <  1)  OR  (X  <=  0.0)  THEN 
BEGIN 
ERROR ; 

EXIT  (CHISQ)  ; 

END 

ELSE  IF  DF  >  40  THEN 
BEGIN 

Y  :=  f  (XPN(X/DF,1. 0/3.0) -1.0)  *2. 0/(9. 0*DF>  I  /SQRT(2.0) 


Y  :=  ( (XPN ( X/DF ,1.0/ 

Y  :=  Y*3. 0* SORT  (DF)  ; 
IF  Y  <  -4.3  THEN 


CHISQ  :  =  1.0 
ELSE  IP  I  >  4.3  THEN 
CHISQ  : *  0.0 

ELSE  CHISQ  :=  Z  (Y,0 . 0, 1. 0)  ; 

END 

ELSE 

BEGIN 

DEN  :*  1.0; 

TEMP  :=  DF; 

REPEAT  .  •  -  -.••• 

DEN  :=  DEN  *  TEMP; 

TEMP  :=  TEMP  -  2.0; 

UNTIL  TEMP  <  2.0; 

POWER  ;*  (DF  ♦  if  DIY  2; 

MUM  :=  XPN(X,  POWER)  *EXP(-X/2.  0)  /DEN; 
IF  ODD  (DF)  THEN 

J  :=  SQRT  (2. 0/X/3.  1415926) 

ELSE 

J  :*  1.0; 

L  : =  1.0; 

M  :=  1.0; 

REPEAT 


DF  :=  DF  ♦  2; 

M  ;=  M*X/DP ; 

L  :=  L+M : 

UNTIL  M  <  5.0000001; 
L  :=  L-  .; 

CHISQ  ;=  J*NUH*L; 
END; 

END;  (*  CftlSQ  *) 
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FUNCTION  INVCHI 


******* ************* **************  ****** 
*  * 

'*  *' 

*  CHI-SQUARE  QUANTILES  * 

*  CASH  * 

*  ALGORITHM  451  * 

*  * 

'*  * 

**************************************** 


CONST 

Cl 

’  = 

1. 565326E-3: 

22 

— 

1 . 06  04  38  E-  3  ; 

C3 

— 

-6.959  356E-3 

C  4 

S 

-1. 32329 3E-2; 

25 

2.2775793-2; 

C6 

3 

-3. 986007E-3 

C7 

s 

-1.51 390 E-2; 

23 

= 

2 . 5300  IE-3 ; 

C9 

= 

-1.45011 7E- 3 

CIO 

SB 

5.  169654E- 3; 

Oil 

3 

-1.15375  IE-2  ; 

C12 

= 

1.  126186  E-2; 

C13 

s 

2.607083E-2; 

014 

S 

-2 .237  35  8E- 1  ; 

C15 

s 

9 . 780499E-5 ; 
-  8 . 553069E-3 

C16 

s 

-8 . 42681 2E-4 ; 

Cl  7 

s 

3.  1255  80E-3 ; 

Cl  8 

s 

C 19 

SB 

1.348028E-4; 

020 

SS 

4. 713941 E-1 : 

- 1 .425296E-2 ; 

C21 

= 

1.0000886; 

A1 

3 

1 . 264461 B-2: 

A2 

s 

A3 

s 

1 .400483E-2: 
-2 . 304527E-2 

A4 

3 

-5.886090E-3; 

A5 

s 

-  1.091 12  IB-2 ; 

AS 

= 

A7 

s 

3.  135411B-3; 

A  8 

s 

-2.728434E-4 ; 

A9 

= 

-9.699681E-3 

A10 

= 

1. 316872E-2; 

A 1  1 

s 

2.6189  IE-2; 
3.483789E-5: 
-8.72971 3E-3 ; 

A12 

= 

-2. 222222E-1 

A 13 

s 

5. 40667E-5 ; 

A  1  4 

s 

A15 

5 

-7 . 27476  IB-4 

A 16 

s 

3 . 2921 81 E-3 ; 

A 1  7 

s 

A 1 8 

S 

4.7140453-1; 

A 19 

= 

1.0; 

VAR 

F,F1,F2,TEMP  :  REAL; 

BEGIN 

IF  (DF  <  1)  OR  (P  <=  0.0) 
BEGIN 
ERROR; 

EXIT  (INVCHI)  ; 

END 

ELSE 

BEGIN 

P  :=  1.0  -  P; 

IF  DF  *  1  THEN 
BEGIN 

TEMP  :*  INVZ  (0.  5*P)  ; 
INVCHI  :*  TEMP*TEMP; 
END 

ELSE  IF  DF  =  2  THEN 
INVCHI  :=  -2. 0*LN (P) 
ELSE 
BEGIN 

?  ; *  DF; 

FI  ;»  1.0/F: 

TEMP  :=  I NVZ  (  1 . 3-P); 
F2  :=  SQRTjtFi)  *TEMP ; 


THEN 


IF  DF  < 
TEMP 


ELSE 

IBM? 


INVCHI 

END; 


(4. 0*ABS  (TEMP)  )  )  THEN 

*F2  +  C2)  *F2+C3)  *?2+C4)  *F2 

♦C6)  *F2*C7)  *PH((  ( ((  (08*C9*F2)  *F2 


mn.o-p)  ; 

TjtFi )  *TEMP ; 

2  ♦  TRONC  (4.0*ABS  (TEMP 
( <<iJ  ( <w1»F2>C2l  *F2*C3 

♦  C5) *F2*C6) *F2*C7)  *? 

♦  CIO)  *F2+C1  1)  *F2*C12 

♦  <JIJ)(*F2*C20)  *P^cil 

(  ((A1>A2*F2)  *P1M  ((A3* 

♦  A5) *?2+A5)  ) *Ff > { (l l 

♦  A10)  *P2*Al  1)  *F2*Al2 


♦  C5)  *F2 *C6)  *F2*C7)  *P1*(  f(U  f08*C9*F2)  *F2 

♦  CIO)  *F2+C1  li  *F2*C12)  *F2*CT3)  «F2*C14)  )  *F1  + 


♦  C 1 7)  *F2*0 1 8)  *F2 


A3*  A4*F2) 


♦  A 1 4)  *F  2*A  1 5)  *F 

♦  A1 3)  *F2*A19: 
TEMPTSM  P*TEMP*F; 


)  >  *? »♦(  (11  (A7*A8*F2)  *P2*A9)  *1 
1  1)  *F 2*  A  i  2)  )  *F1*  ( ( (  (  (A13*F2 
15)  *F2*A16)  *P2*Ai7)  *F2*F2 
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FUNCTION  PI 


************************  **************** 

*  * 

*  * 

*  P  DISTRIBUTION  *! 

*  SMALL  DEGRESS  OP  FREEDOM  * 

*  STATISTICAL  COMPUTING  * 

*  P.  114  * 

*  * 

*  * 
**************  ********************  ****** 

CONST  PI  =  3.  14  159265; 

VAR  A. A 1. Ail. AN2, ART, D1 ,  D 2 , D3 ,R ,S 1 , S2 , DEL ,XM, XK, T , 
TEM,C  ;  REAL: 

I#M,N  ;  INTEGER; 

BEGIN 

AN1  :=  DPI; 

AN 2  :=  DP2; 

A  :=  DP 1*X/ (DPI *X+DF2) ; 

A 1  :  =  1  -  A: 

IP  A1  <=  0.0  THEN 
A 1  :=  1. OB-37; 

D1  :=  AN1*0.5 ; 

D2  :=  AN2*0.5; 

D3  ; *  D 1  ♦  D2  -  1.0; 

R  : =  0.0: 

51  :=  0.0; 

52  :*  0.0: 

DEL  :=  1.5; 

XM  :*  1.0; 

XK  :=  1.0; 

C  : *  0.25; 

N  :=  DP2; 

NHILE  C  <  0.  875  DO 
BEGIN 

M  :=  TRUNC  (  D2)  ; 

M  :=  2  *  M; 

IP  a  -  N  THEN 
BEGIN 

N  :=  TRUNC  (02)  -1 ; 

IP  N  >  0  THEN 
BEGIN 

FOR  I  :=  1  TO  N  DO 
BEGIN 

51  :=  DEL*S1*R; 

D2  :=  D2  -  1.0; 

D3  :=  D3  -  1.0; 

TEM  :*  A1/D2; 

R  :=  D3*TEM; 

52  :=  (R+TEMMS2; 

END; 

END; 

51  :*  DEL+S  1*R; 

DSL  :■  0.0; 

T  :*  -1.0 : 

D3  :  *  —1.0; 

52  :*  A*S2; 

C  :  *  C  ♦  0.5; 

END 

ELSE  (*  DEGREES  0?  FREEDOM  ODD  *1 
BEGIN 

N  :*  TR  UNC  (D2)  ; 

IF  N  <>  0  THEN 


jMujlilTI  l»  "  -  .•  “ 
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FUNCTION  P2; 


******************  ***************************** 

*  P  DISTRIBUTION  * 

*  LARGE  DEGREES  OF  FREEDOM  * 

*  EDUCATIONAL  AND  PSYCHOLOGICAL  MEASUREMENT  * 

*  7.  25,  SO.  3,  P.  877-879  * 

**************************************  ********* 


CONST 

Cl  *  0. 19685'4  ;C2  * 
7AR 

S,T,Z, Y,XX,TY,J,K  i 


0. 115194;C3 
REAL; 


0.000344;C4  =  C 


BEGIN 
IF  X  > 

BEGIN 

s  :  = 

T  :  = 

Z  :  = 

END 
ELSE 
BEGIN 
S  :* 

T  :  = 

Z  :  = 

END; 

J  :»  2.0/9.0/S; 

K  •  *  2  0/9  0/T* 

T  ;*  ABS7(i.0-ft)*XPN(Z,  (1.0/3. 0)  »  -1.0+J)/ 


1  THEN 

DF1 ; 
DF2; 

x; 


DF2; 
DF1 ; 
1.0/X; 


T  ;*  ABS((1.0-K)*XPN(Z,  (1.0/3. 0H-1.0*. 
SQRT(K*XPN(Z,  (2.0/3.  j[[  *J)  ; 

if  ;*  i.o  ♦  r*(ci*r*(C2+7*(C3*i*c4n )  ; 

XX  :*  0.5/XPN  (YY,4.0|  ; 

IF  X  >»  1.0  THEN 
F2  :»  1.0  -  XX 
ELSE 

F2  :*  IX; 

END;  (*  F2  *) 

FUNCTION  F; 

BEGIN 

IF  (DPI  <  1)  OR  (DF2  <  1)  OR  (X  <  0.0)  THEN 
BEGIN 
ERROR; 

EXIT  (F)  ; 

END; 

IF  ((DPI  <  100)  AND  (DF2  <  100))  OR 
(DPI  <  20)  OR  (DF2  <  20)  THEN 
'  F  J*  FI  (X , DF  1 , DF2) 

ELSE 

F  :*  F2  (X , DF 1, DF2) ; 

END;  (*  F  *) 
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FUNCTION  INVF  ; 

, *************************************** 
* 

*  P  OUANTILES 

*  BISECTION  SEARCH  FOR 
SHALL  DEGREES  OF  PRSEDOM 

*  LARGE  DEGREES  OF  FREEDOM 

*  HANDBOOK  OF  MATHEHAPICAL  FUNCTIONS 

*  P.  9t»7 


******** 


******************************* 


CONST  EPS  *  0.005; 

*A?EHP, PT, ENDR, ENDL, «IDPT,STSP,N,H,T,r,Z  :  REAL; 

BEIFS(DF1  <  D  OR  (DF2  <  1)  OR  (P  <=  0.0)  OR  (P  >=  1. 
BEGIN 
ERROR; 

EXIT  (INVF)  ; 

IlSE  IF  (DF1  =  1)  OR  (DF2  =  1)  THEN 
BEGIN 

IFINV?  J-1s5H!ifTn0-3.(1.0.W  >  * DF2>  ) 

SLllvF  :  =  1.0/SQR(INVT((1.0-P/2.0)  ,DP1))  ; 

END 

ELSE 

BEzI:»  INVZ(P): 

I  M  J?0Z*‘i:5/)(<?5/;(DF2;1.0)  ♦  .hO/CDrirl.OIL; 

"  inal:i/6:V  -2.o/H:o^?f!*  ”  ' 

IF:(DF1P^=*  15)°  da  (DF2  <=  15)  THEN 
BEGIN  „ 

STEP  :=  0 . 1*T; 

ENDR  ;*  T; 

REPEAT 

TEMP  :=  ENDR; 

ENDL  :=  ENDR  -  STEP; 

IF  ENDL  <=  0,0  THEN 
ENDL  :=  0.0; 

ENDR  *  s  ENDL4 

UNTIL^PT^i3??^^}^!^*-  ENDL)  <=  BPS); 

IFDAbI(3NDR?-  ENDL)  >  EPS  THEN 

BE»HILS  (ENDR  -  ENDL)  >0.1  DO 

BE?r:*  (ENDR  ♦  ENDL) /2. 0; 

PT  :*  FI  (T.DF1  ,DF2)  ; 

IF  PT  <  P  THEN 
ENDL  :*  T 
ELSE 

ENDR  :*  T; 

END; 

END; 


0)  THEN 


0))* 


APPENDIX  B 


PROCEDURE  CONVERTSTRIN3 ; 

*************************************** 

*  * 

*  PROCEDURE  TO  CONVERT  A  STRIN3  * 

*  VARIABLE  TO  A  NUMBER  * 

'*  * 

************************  *************** 

LABEL  1; 

VAR 

PLACEINT.SIGN.I  :  INTEGER; 

BFDEC  :  BOOLEAN; 

ACC, TACC, PLAC E  :  REAL; 

PREVCH, CH  :  CHAR; 

BEDATA  :*  CONCAT  (DATA,'  ')  ; 

prevch  ; 

BFDEC  ;»  TRUE; 

INT  :3  0: 

RL  :*  0.0: 

ACC  :=  0.0; 

PLACE  ;=1; 

FORHi:T*  1  TO  LENGTH  (DATA)  DO 
BEGIN  _ 

CB  :=  DATA  I  ; 

IF  CH  =  PREVCH  THEN 
GOTO  1; 

l|2CHHirC'^i.'9','.',*  'r*3  THEN 

CASE  ^  CH  ^OFf  j  ^  j  j  -  j  |^|  |  <j  |  ,  g  |  i^v. 

BEGIN  * 

IF  BFDEC  THEN 

BERLS:=  RL  *  PLACE  ♦  (3RD  (CH)-48)  ; 

PLACE  :=  10; 

END 

ELSE 

BEGIN 

llk?.  il  SL|SScS  *-li.D(CW-«8, 

END; 

END; 

BEGIN 

BFDEC  :=  FALSE; 

PLACE  : =  1 

END; 

«  •  • 

BEGIN ^ 

IF  IflS^RL)  3< ^3^767 .0  THEN  <*  HAX  IN 
INT  :*  TRUNC  (RL)  (*  IS  327 


ELSE 

INT 


i; 


BEGIN 

SIGN  :*  -Is 
PRE7CH  :=  ' 
END: 

'  ' :GOIO  1; 

END*  (*C*SE*) 

END: 

(*  C0N7E RISTHI NG 


*) 


na 


LIST  OP  BEPEEBEHCES 


1. 


2. 


3. 


4. 

5. 
S. 
7. 


iaa  af 


ft 


asters 


Q.L.QGQA- 
es  is. 


Ling,  a.  F. ,  General  Considerations  on  the  Design  of 
an  Interactive  Systam  for  Data  Analysis,  Cqeean..i ca¬ 
tions  of  the  ACM,  v.  23,  p.  147-154,  3 arcfl  l 9 oOT 


Bullae,  B.  £• #  • 
Packages  for  the 
Statistician.  ?. 


1§30  *  s 
34 ,  no 


of  Statistical  Computing:  ihat 


iSilS:  l: 


Kellner,  J.,  “Pascal  Operand  Formats,1'  Apple  Orchard, 
v.  1,  no.  2,  p.  39,  1930. 

Forsythe,  A.  B. ,  “Elements  of  Statistical  Computation," 
Byte,  p.  182-184,  January  1979. 


Dixon,  B.  J.  and  Bassey,  F. 

8d-'  ”• 


3.  Kennedy ,  2.  J.  and  Gentle,  J.  E.,  Statistical  Compu¬ 
ting.  p.  112-116,  Dekkar,  1980. 

9 . ,  Johnson,  H.  L.  and  KotZ,  S.x  Continuous  Univariate  Pis 
tributions.  v.  2,  illey,  1970. 


10. 


Storer,  8., 
tor .  aasters 


11. 


A|^le  Computer,  Inc., 


Apple  Pascal  Reference  Manual. 


12. 


Gibbons, 
tion  and 
v.  29,  no 


J.  and  Pratt,  J.  ».,  “P- Values:  Interpreta- 
Hethodology^  T^e^Americ^^|tatistician. 


49 


BIBLIOGPAPHI 


Abraaowitz,  M. 
Sanctions,  p. 


and  Stegun,  I.  A.,  Handbook  afc  jt a 
127-949,  Do7er  Publications, 19/9. 


atheistical 


Association  of  Computing  Machinery,  “collected  Algorithas 
froa  ACM,  1980. 


Dixon,  i.  J.  and  Massey,  P.  J.,  Introduction  ta  Statig- 
^  ' — , — ■=  -  125-140,  Osborae/McGraw-Hill, 


tical 
Inc. , 


sis,  3rd  ei. ,  p, 


Golden,  B.  B. ,  Weiss.  D.  J.,  and  Dawiss,  a.  7. ,  ‘An 
Bvaluation  of  Jaspen*  s  Approx  jjaation  of  the  Distribution 
Functions  of  the  ?,  T,  ana  Chi-sguare  statistics, 
l^ca^ig^al  and  Psychological  Measurement.  7.  28,  p.  163- 


Johnson,  M.  L.  and  Kotz,  S..  Continuous  Uni7a.riate  Pis 
trihut ions.  7.  1,  Wiley,  1970. 


Mood,  A,  M,  and  Graybill,  F.  A...  Introduct 
of  statistics.  2nd  ed. ,  McGraw-Hill,  1963. 


ion  to  thP  Theory 


Poole,  L.  and  Borchers,  M.,  Sale  goaaoa -BASIC.  Proqiaas,  3rd 
ed.,  p.  125-140,  Osborne/McGraw-Hili; 

Saillie,  K.  W.  and.Anstey,  T.  H.."A  Mote. on  the  Calcu- 
*  ion  of  probabilities  in  an  F-Distribution,  Canauni- 

iona  of  the  Association  of  Coaoutina  MachiaarY,  ▼  .  7,  p. 


Zehna,  P.  W.,  introductory  statistics,  p.  191-280,  Prindle, 
Weber  and  scnaidt ,  Inc. ,  1974 . 


i 


t 


initial  distribution  list 


2. 


3. 


4. 


5. 


6. 


Defense  Technical  Inforaation  renter 
Caaeron  Station 
Alexandria,  Virginia  22314 

Library,  Code  0142 

Naval  Postgraduate  School 

Monterey,  California  93940 

Professor  D.  H.  Baer.  Code  55Bn 
Department  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  California  93940 

LCDR  C.  F.  Taylor,  Code  55Ta 
Departaant  of  Operations  Research 
Naval  Postgraduate  School 
Monterey,  California  93940 

MAJ  R.  a.  Duff 
9  Revere  Road 

Monterey,  California  93940 


MAJ  j.  D.  Morgeson 
903  Virginia 
Crane,  Texas  79731 


No.  copi 
2 


2 


1 


1 


1 


2 


i 


i 


5  1 


♦ 


